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Spectral matching for abundances and clustering analysis 
of stars on the giant branches of uj Centauri 
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ABSTRACT 

We have determined stellar parameters and abundances for 221 giant branch stars in 
the globular cluster u) Centauri. A combination of photometry and lower-resolution 
spectroscopy was used to determine temperature, gravity, metallicity, [C/Fe], [N/Fe] 
and [Ba/Fe]. These abundances agree well with those found by previous researchers 
and expand the analysed sample of the cluster, fc-means clustering analysis was used 
to group the stars into four homogeneous groups based upon these abundances. 

These stars show the expected anticorrelation in [C/Fe] to [N/Fe]. We investi- 
gated the distribution of CN- weak/strong stars on the colour-magnitude diagram. 
Asymptotic giant branch stars, which were selected from their position on the colour- 
magnitude diagram, were almost all CN-weak. This is in contrast to the red giant 
branch where a large minority were CN-strong. The results were also compared with 
cluster formation and evolution models. Overall, this study shows that statistically sig- 
nificant elemental and evolutionary conclusions can be obtained from lower resolution 
spectroscopy. 
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1 INTRODUCTION 

to Centauri is a large star system in the southern sky. Its 
overall appearance is that of a globular cluster (GC) but 
it has several features that set it apart from the classi- 
cal "simple-stellar population" globular cluster. Its main 
sequence is at least bifurcated ( |Bedin et al.| |2004| |Pi-| 
|otto et aL 2005 I and the giant branch (GB) shows multi- 



ple branches (Johnson & Pilachowski 2010 and references 



— 1.50, with evidence for further peaks at —1.15, —1.05 and 
—0.75 (J&P10). These peaks correspond to different GB 
that have been observed photometrically ( |Rey et aL||2004[ ). 

uj Cen differs from most mono-metallic clusters not only 
in its metallicity range but in its neutron-capture abun- 
dances. MoslQ mono-metallic clusters show no variation in 
their s- or r-process abundances (e.g. in 47 Tuc as observed 
by Worley & Cottrell||2012 1, while within oj Cen there is a 



therein) (herein referred to as J&P10). These multiple path- 
ways on the colour-magnitude diagram are chemically dif- 
ferent in overall metallicity (Norris & Da Costa 19951 and 



helium abundances ( jDupree et al.|2011[ ). |Marino et al.| ( |2012[ ) 

(herein referred to as M+2012) suggest that some of these 
difference could be explained by CNO variations. 

A number of studies have examined the metallicity and 
abundance profile of uj Cen in detail. Cluster stars with 
metallicities from [Fe/H]= —2.2 to —0.5 have been observed 
spectroscopically. Most of the GB stars are found to be 



positive correlation of s-process abundances with metallic- 
ity. In fact there is over a 1 dex increase in abundance of 
[La/Fe] and [Ba/Fe] between [Fe/H]= -1.9 and -1.5. Bar- 
ium abundances from | Villanova et al.|p010 l have a 1.5 dex 
scatter at all [Fe/H]. This was also found for the light s- 
process elements, yttrium and zirconium. |Villanova et al.| 
( |2010 1 reported evidence for bimodality in the [Ba/Fe] to 
[Fe/H] distribution. For stars with [Fe/H]< —1.5 there were 
two populations: [Ba/Fe]-~ +1.0 or ~ —0.2. This result does 
not agree with the larger samples of |JfcP10| and |Marino| 



around [Fe/H]= -1.75 


Norris & Da Costa||1995 Stanford 


et al.|2006||Marino et al. 


201 1 Johnson & Pilachowski 2010 1 . 



et al. (20111 who found no such bimodality. Instead there 



was a continuous trend of increasing [Ba/Fc] with increasing 



There is a second peak in the metallicity distribution round 
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1 One mono-metallic cluster with a neutron-capture spread is 
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[Fe/H], at least up to [Fe/H]= —1.4. Above this metallicity 
there is a flattening or even turning over of [Ba/Fe] with 
respect to [Fe/H] ( jStanford et~aL|2010| ). 

There are a number of small studies of the carbon 
and nitrogen abundances in uj Cen. Brown & TWallerstein| 
(1993]) analysed six stars finding -0.45 <[C/Fe]< +0.1 and 
0.75 <[N/Fe]< 1.25. These stars had a large range of metal- 
licities but did not include any stars with [Fe/H]> —1.25. 



Table 1. Overview of different data sources used in this work. 



Norris & Da Costa ( 1995 ) had a larger sample of 40 (biased) 



stars for which they determined a larger range than seen by 



Brown & Wallerstein ( 1993 1 for both carbon and nitrogen. 



They found that there was a "floor" to their [C/Fe] values 
at ~ —0.8, a feature that they pointed out was not due to 
their abundance analysis method. Their most carbon-rich 
objects were roughly solar. jStanford et al.| ( |2010[ ) examined 
33 stars at low resolution again finding a large range of car- 
bon and nitrogen values in the cluster. They reported an 
anticorrelation but no bifurcation in their results. 

The most recent study and at the highest resolution 
(~ 20, 000) is|M+2012| They observed 77 giants and found a 



correlation between the overall C + N + O abundance and the 
metallicity. Using the sodium-oxygen anticorrelation, they 
were able to distinguish what they referred to as first and 
second generation stars. Their results show a clear bifurca- 
tion between CN-weak and CN-strong stars (their Figure 
2). The nitrogen-rich stars are grouped together with little 
spread in carbon. The CN-weak stars show a range of [C/Fe] 
from —0.5 to +0.5. With 221 stars, the work presented in 
our study will expand upon these stars as there are only 10 
stars in common with |M+2012| 

One of the important ingredients in any theory of the 
formation of uj Cen is the age spread of the stars. A va- 
riety of ranges have been reported in the literature. The 
shortest age differences have been reported by|Sollima et al 



( 2005 1 who found 2 Gyr. At the opposite end is Villanova 
|et al.| p007[ ) who found up to 5 Gyr. The placement of the 
isochrones that arc used in this sort of ago determination de- 
pends on the CNO abundance of the models. |M+2012"] used 
the |Pietrinferni et al.|p006] ) isochrones and found that "for a 
fixed turn-off and sub-giant branch (SGB) brightness, CNO- 
enhanced isochrones provide younger ages than isochrones 
corresponding to a canonical a-enhanced mixture; as a rule- 
of-thumb we found the <5age/<5[CNO] ~ — 3.3Gyr per dex". 

Multi-object spectrographs allow for hundreds of stars 
to be observed in one night, making it possible to efficiently 
sample an object like a globular cluster with 10 4 — 10 6 stars. 



Recent papers such as |Marino et al. (|2011| ) used the VLT 
to obtain spectra of 300 stars and |JfcP10| published results 
for 855 stars. They confirmed the peculiar chemical make- 
up of uj Cen. Our study makes use of both of these datasets 
and extends them with an investigation of the |van Loon] 
et al. (20071 (herein referred to as vL+07) spectral library. 
The work presented here determined the common stellar pa- 
rameters (temperature, surface gravity, metallicity) as well 
as s-process, carbon and nitrogen abundances using lower- 
resolution spectra and photometric data of this sample of 
uj Cen giant branch stars. 

Section [2] discusses the observational data sources; Sec- 
tion [3] explains the methods for analysing the |vL+07| spec- 
tral library; Section [4] presents the results and quality as- 
surance; Section [5] presents a statistical clustering analysis 
of the dataset from this work and of Marino et al. (2011 



Source 



Type 



Stars 



vL+07 



J&P10 



Bellini ct al. 

2MASS 



(2009) 



M+2012 



Spectroscopy 

Abundances 

Photometry 

Photometry 

Abundances 



1518 
855 
359,391 
26,985 
77 



2012{ ); Section [6] shows that the easily identifiable asymp- 
totic giant branch (AGB) stars are likely to be CN-weak; 
Section[7]discusses the results in the context of evolutionary 
models proposed in previous works. 



2 OBSERVATIONAL DATA SOURCES 

A dataset of 221 giant branch stars (including 14 stars that 
were observed twice spectroscopically) was created by com- 
bining information from |vL+07 2MASS 
[2006] ), 



Bellini et al. (20091 and 



Skrutskie et ah] 
JfcP10| Table [1] describes 



the data sources that were combined and used in this work. 

The datasets, and how they were combined are de- 
scribed in this section. Section |2.1| describes the |vL+07| 
spectral library, Section [2.2| describes the photometric data 
sources, with the positional matching of different datasets 
described in Section [23] The final selection of stars is found 
in Section [2~4l 



2.1 |vL+07| spectral library 



The |vL+07| dataset is a spectral library of 1518 post-main 
sequence stars in uj Cen. Spectra were obtained with the 
2dF instrument at the Anglo-Australian Telescope, cover- 
ing approximately A ~ 3840-4940 A at a resolving power of 
A/AA ~ 1600 and a signal-to-noise per pixel ranging from 
~ 50 in the blue to > 100 in the redder part of the wave- 
length range. For each st ar |vL+07| provided a unique ID 
from |van Leeuwen et al. ( 2000 1 (known as the LEID) , po- 
sitional information, a B magnitude, a B — V colour and 
a radial velocity. |vL+07| also did a stellar parameter de- 
termination for their catalogue. This was done solely from 
the spectra and assumed solar-scaled abundances. For our 
study, we decided to use their raw spectra and analyse the 
stars ourselves with appropriate globular cluster abundances 
and concentrating on small sections of the spectrum. 

We iteratively continuum normalized the raw spectrum 
by mapping the spectrum onto a synthetic spectrum of the 
best stellar parameters. Three anchor points were selected: 



4088A, 42 20A and 4318A (as also used in |Worley fc Cot- 
trell||2012 1. These three points were joined by two straight 



lines and the intensities were mapped from the raw spec- 
trum onto the synthetic spectrum. This method is depen- 
dent on the synthetic spectrum being used and the contin- 
uum placement was recomputed at each iterative step using 
the method described in Section \S. 21 

At the resolution of the |vL+07| spectra, and with 
its wavelength coverage, it was not possible to use a 
spectral-temperature indicator such as H/3. Therefore the 



temperature-colour relationships of Alonso et al. ( 1999 1 were 
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used. This could then be combined with the luminosity re- 
lationship to estimate the surface gravity of the stars. For 
stars of unknown metallicities, the most appropriate colour 
to use is V—K, as noted by Alonso et al. ( 1999 I in their sec- 
tion 2.3. Consequently, additional photometry of the stars 
in the |vL+07| catalogue was needed, since they only had B 
and V. 



2.2 Photometry data sources 

2.2.1 2MASS photometry 



The Two Micron All Sky Survey (2mass) (Skrutskie et al 



20061 observed the sky in the near-infrared covering three 
photometric bands: J (1.25 /im), H (1.65 fim), and K s 
(2.17 fim). The 2MASS data were sampled for a circular 
region with a radius of 28" around RA = 13 h 2 6 m 14.16 5 , 
Dec = -47°31'11.95", the midpoint of the |vL+07| data. This 
created a catalogue of 26,985 stars with J, H and Ks mag- 
nitudes in the 2MASS photometric system. The majority of 
the stars in the raw dataset were from the globular cluster 
with obvious contamination by field stars. The photometry 
extended well below the cluster's main sequence turn off, 
covering all the [vL+07] stars. 



2.2.2 



Bellini et al. \2009 ) photometry 

Bellini et al. (2009) used ESO archive data to create a 
new, CCD-based, proper-motion catalogue for ui Cen, ex- 
tending to B ~ 20, containing 359,391 stars with U, B, V, 
Rc, Ic, and Ha. They used 279 archive images acquired at 
the ESO/MPI2.2m telescope at La Silla (Chile) equipped 
with the wide-field imager camera (WFI). It covered 83% 
of the stars of the |vL+07| dataset. This catalogue improved 
upon the previous large proper motion catalogue of uj Cen, 



produced by van Leeuwen et al. (2000) from photographic 



plates. Although |vL+07| had B and V magnitudes for all of 
their stars, it was decided to use the |Bellini et al.| ( |2009[ ) 
photometry instead as it provided a very precise colour- 
magnitude diagram that could be used to distinguish the 
metal-poor AGB (the bluest stars on the GBs)from all of the 
red giant branch (RGB). Other AGBs will be mixed with the 
RGBs as a consequence of the different evolutionary paths 
followed by stars of different ages and metallicities. 



2.3 Positional matching 

To determine the T e ff and log g of the stars, the V—K colour 
was needed. This required the positional matching of vL+07| 



(2009) 



sources and the photometric sources of Bellini et al. 
and 2MASS. The closest positional match was found. 

2MASS stated that their positional accuracy should be 
less than 130 mas for saturated sources and < 80 mas for 
the best unsaturated sources. The van Leeuwen et al. (2000) 
dataset, from which the |vL+07| is derived, had positional 
errors mostly less than 100 mas in right ascension and dec- 
lination. Bellini et al. ( 2009 1 stated that based upon ~ 5500 
reference stars, their error was 45-50 mas in each coordinate. 
The positional matching of the three photometric catalogues 
was done using TOPCAT ( |Taylor|2005[ ). 982 stars were found 
to be in common between |vL+07| and 2MASS, with 672 of 



Table 2. Photometry of the 221 stars. The LEID is the star 
label given by |vL+07| The U and V photometry is from [Bellini| 
|et al.||2u09| and the J and K photometry is from 2 MASS. All the 
photometry is in their original system. The number of spectra of 
each star in the |vL-)-07| spectral library is given. 



LEID 


Tja 




J 6 


K b 


N. of Spectra 


15022 


14.19 


12.10 


9.730 


8.95 


1 


15023 


14.15 


12.26 


9.964 


9.23 


1 


16019 


14.78 


13.28 


11.268 


10.62 


2 


16027 


14.83 


13.55 


11.532 


10.86 


1 


17029 


14.42 


13.17 


11.201 


10.54 


1 


17046 


14.41 


13.31 


11.364 


10.74 


1 


18040 


15.10 


12.79 


10.448 


9.60 


2 


19022 


14.30 


13.16 


11.212 


10.58 


2 


19062 




12.85 


10.601 


9.87 


1 



(This table is available in its entirety in a machine-readable 
form in the online journal. A portion is shown here for guidance 
regarding its form and content.) 




Figure 1. A V-(U — V) colour-magnitude diagram of u Cen de- 
rived from the |Bellini et al.| l |2009| photometry (small dots) . The 
open squares highlight the positions of 172 (of 211; the missing 
stars did not have U photometry) stars in our sample. The magn- 
tiude limit of V ~ 13.5 is due to the limiting magnitude of |JfeP10| 
The abundances determined are presented in Table [4] 



2.4 Selection of stars 



|JfcP10| determined several elemental abundances including 
[Fe/H], [O/Fe], [Na/Fe] and [Ca/Fe]. For the analysis pre- 
sented in this work it was decided to concentrate on stars 
that were also in the |JfcP10| catalogue as it would result in 
more precise abundances that did not have further uncer- 
tainty due to unknown oxygen and calcium abundances. A 
full analysis of the complete |vL+07| spectral library, with 
reasonable estimates of oxygen and calcium, is in prepara- 
tion. 

IJfcPlOl selected their stars from [van Leeu wen et al.l 
pOOO]) and there were 291 matches in total between |JfcP10| 
and |vL-|-07| Not all of these stars were positionally matched 



to the 



Belli ni et al.| ( |2009[ ) and 2MASS libraries. Limiting to 



these stars also matching to Bellini et al. (20091. 



just those stars that were, a final catalogue of 221 stars was 
created (Figure [l]and Table [2). 

This dataset does have a selection bias against the core 
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of the cluster. This is due to the crowding that makes it diffi- 
cult for a best positional match to identify the correct star. 
Comparing the radial distributions of the J&P10 dataset 
and our own, we find that their stars are much more cen- 
trally concentrated: 57% of the stars in the |JfeP10| dataset 
are within 6' of the cluster centre while for our set it is only 
29%. In Section |5.3| we note that this selection effect does 
not appear to have greatly biased our sample, as we have a 
similar proportion of O-rich and O-poor stars as found by 
IJfcPlOl 

Overall, our 221 stars have a magnitude range of 11.2 < 
V < 13.6 and a colour range of 0.8 < (B — V) < 1.8. This 
corresponds to stars with T e s from 3800 K to 4900 K and 



logg from 0.3 to 1.8 (Section 3.1 1 



3 DATA REDUCTION OF SPECTRAL 
LIBRARY 

3.1 Temperature and gravity 



The temperature-colour relationship of Alonso et al.| ( |1999[ ) 
was used to calculate the T e a using their equations 8 & 9. 
These equations have constraints on V — K and [Fe/H]. If 
the star satisfied the constraints for both equations then an 
average of the two temperatures was used. If neither set of 
constraints were satisfied then the star was excluded. 

Alonso et al. (1999) required the V — K photometry to 



be in the Carlos Sanchez Telescope photometric system. The 



photometry of Bellini et al. (2009) and 2MASS photometry 



was converted to this system using a method developed by 



Johnson et al. (20051. The interstellar reddening adopted 
was E {V-K) = 0.3 \ Johnson et al.||2008[ ). |Calamida et aL 



( 2005) found evidence for an intra-cluster reddening varia- 
tion in u Cen of 0.06 sC E(B-V) ^ 0.13. In terms of the 
effect on the T e g, this would create a range in T a g of ~ 70K 
for the hottest stars. Such a change in temperature is of a 
comparable order of magnitude to the uncertainties result- 
ing from the photometric errors. 

Differential reddening was disputed by |Villanova et aL] 
(2007[). They produced CMDs of a million stars from ACS 



photometry. They point to the sharpness of their CMDs as 
evidence for the lack of any range of intra-cluster reddening 
for to Cen. So although the range of reddening suggested by 
Calamida et al.| ( |2005| l would have a measurable effect on 



the temperatures, it did not seem to be a significant issue 
across the whole cluster. 

The surface gravity was determined using the elemen- 
tary equation, 



log g, 



with M, 



O.4O(M bo i > .-M b ^0+Iogfl0 
+41og(T»/T ) +log(M./Af ), 



+4.75, T = 5777 K and log fl0 = 4.44 



(Allen & Cox 20001. The stellar mass was set to 0.8 M . 
There is evidence for an age range in uj Cen which would 
imply a mass range. This mass range would extend down to 
0.6 Mq. However, since the logarithm of the stellar mass 
is used here, this reduced the effect on the final gravity 
value to 0.1 dex at the most extreme. As such, the prob- 
able mass range was ignored. The bolometric correction was 



dependence which, across the metallicity range of uj Cen 
would affect the surface gravity by 0.01 dex. 

Differences in gravity can be mimicked by differences in 
metallicity, especially with the ionized barium line at 4554A 
used in this study for determining [Ba/Fe]. Increasing the 
gravity of a star will decrease the strength of spectral lines 
of ionized species (for all other parameters remaining con- 
stant). It was found that a 0.3 dex change in gravity (for a 
T cff = 4250K, log g = 1.6 and [Fe/H]= -0.5 star) could be 
compensated for by a 0.05 dex change in the metallicity. For 
a gross change of metallicity it was not possible to use the 
gravity to compensate. 

| JfcP10| and |Marino et "al] |20TT| [2012] ) all found essen- 
tially the same temperatures and gravities for our sample of 



stars as they used the same method ( Johnson et al. 2005 \ 
and photometry ( |van Leeuwen et al.||2000| |Bellini et al. 
2009 1 . Any small differences were due to distance moduli, 



reddening and stellar mass selections. 



3.2 Spectral matching 



Spectrum synthesis was performed using MOOG (Sneden 



19731. MOOG is a code that performs a variety of lte line 



analysis and spectrum synthesis tasks. Around this was cre- 
ated a Python wrapper to find the best matching synthetic 
spectrum using known parameters and physically reasonable 
estimates for unknown parameters of the stars. 

For each star, the temperature and gravity were calcu- 
lated from the V — K photometry and assumptions about 
the stellar mass as described in Section 13.11 The metallic- 
ity was also known for all of these stars from |JfcP10| To 
keep the results of [Fe/H], [C/Fe] and [N/Fe] self-consistent, 
[Fe/H] was determined for all the stars using our spectral 
matching method. Being able to determine [Fe/H] will be 
required when working with the full |vL+07| dataset as the 
metallicities are not known for the majority of these stars. 

At the resolution of the |vL+07| spectra, there was not 
enough information to try and determine an [O/Fe] abun- 
dance, so the |JfcP10| values were used in this work. The 
other known value for all of the stars was [Ca/Fe] which 
were also taken from lJfePlOl 

The T e s, log g and [Fe/H] values were used to generate 
the necessary model atmospheres via interpolation 2 ] within 
the a-enhanced grid of Castelli & Kurucz ( 2004 1. A value of 



[a/Fe]= 0.3 was used ( |JfcP10[ ). For t he CH and CN region 
of the spectrum, line lists were from Norris (20121. These 



line lists were found to match to the high resolution atlas of 
Arcturus ( Hinkle et al.|2000 i using abundances determined 
by Decin et al. (2004). They also matched the Arcturus spec- 



trum convolved to the resolution of the | vL+07| spectra. The 
atomic log<?/ values in the line list for the barium region 
investigated (400 to 500 nm) were adjusted so that the syn- 
thetic spectrum matched the Arcturus spectrum, again at 
both resolutions. 

Lines of barium exhibit hyperfine splitting (hfs). In 
spectrum synthesis, hfs will cause the line to have a larger 
equivalent width for a given abundance compared to when 
hfs is not included. For the 4554A line of barium, two line 



taken from Alonso et al. (19991 and there is a small [Fe/H] 



2 The interpolation routine 

https: / / github.com /andycasey / atmosphy. 



used 



atmosphy 
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Figure 2. Examples of three spectra, (a) is an RGB-a star with 
a high metallicity ([Fc/H]= —1.0), illustrating the three regions 
that were used in to determine the metallicity (via the Ca 4224A 
line), carbon (CH) and nitrogen (CN) abundances. The noise level 
is shown to reflect the signal-to-noise in the spectra, (b) and (c) 
are a CN-weak/-strong pair and show the variation in the CN 
band. The difference in [N/Fe] between (b) and (c) is 1.4 dex. 



lists were created: one that included all the isotopic species 
of the barium feature (Ruttcn 1978 1 and one that did not. It 
was found that across the temperature and gravity ranges of 
the cluster, the effect was negligible. Such results confirmed 



what was found by Stanford et al. (2010) for main sequence 



stars. They found a difference of [Ba/Fe]< 0.05 dex between 
abundances found with hfs and those without. 

The microturbulence (vt) was assumed to be 2.0 kms -1 
for all of the stars. This is within the range of values deter- 
mined by |JfcP10| We were unable to refine it further due to 
the lack of unblended iron lines to remove the trend in iron 
abundance with line strength. 

The code we developed read in the spectrum of the star 
along with previously determined stellar parameters (T e g, 
logy, [O/Fe] and [Ca/Fe]) for the star and then found the 
best fitting synthetic spectrum created using MOOG. It was 
designed to search the parameter space in a systematic way. 



The S(3839) index from vL+07 was used as an estimate as 
to whether the star is CN-weak and CN-strong. This gave 
the initial guesses for [C/Fe] and [N/Fe]. An iterative pro- 



cess was used. At each step of the process a model atmo- 
sphere and then synthetic spectrum was created from the 
full parameter set (T cff , logy, [Fe/H], [0,Ca,C,N,a/Fe]). The 
[Fe/H] was determined using the metallicity sensitive Cal 
line at 4226A. The other two spectral features used were the 
carbon-dependent G band (mainly CH: A 2 A — X 2 H) region 
(4295A to 4325A) and the CN blue-system (£> 2 E - X 2 S) 
(4195 A to 4222 A). 

At the first step, the initial guess for [C/Fe] and [N/Fe] 
was used to find [Fe/H] using the Cal and the known [Ca/Fe] 
ratio (J&PlOj), minimizing the difference in the equivalent 
width between the synthetic and observed spectra for the 
calcium line. Then the newly determined [Fe/H] was used in 
the x 2 minimization to find [C/Fe] from the G band. Finally 
the determined [Fe/H] and [C/Fe] were used to determine 
[N/Fe] from the CN-region. This process was then repeated, 
with the newly determined parameters from the previous 
iteration. [Fe/H] was redetermined using the new [C/Fe] and 
[N/Fe] values, requiring the raw spectrum to be continuum 
normalized to the new synthetic spectrum. This full cycle 
was repeated until it converged on an unchanging triplet 
of [Fe/H], [C/Fe] and [N/Fe]. A sample of three spectra is 
shown in Figure [2] 

A delta analysis was performed to understand the sensi- 
tivity of the spectral matching method to the input param- 
eters of temperature and gravity. Four stars were selected, 
one from each group identified in Section [5] with two stars 
being low temperature (~ 4000 K) and two being higher 
(~ 4500 K) . The results of this analysis are presented in Ta- 
ble^ It was found that increasing the temperature by 100 K 
results in 0.1-0.2 dex increase in [Fe/H] and [C/Fe]. [N/Fe] 
changed by ±0.1 dex except for the most N-rich star which 
decreased by 0.4 dex. Changing the gravity by +0.3 dex had 
little to no effect on the abundances (see also comments in 
Section 3.1 1. 



The first 20 stars (which have a total of 25 spectra) 
were analyzed by repeating their analysis with four random 
starting positions. This was to test the convergence of the 
method. All but three of the stars returned the values af- 
fected by less than 0.1 dex. The other three stars were all 
very N-poor stars and were among a group of stars that were 
excluded from the comparative analysis. At the spectral res- 
olution of the [vL+07] library, there is not a lot of sensitivity 
to nitrogen at low values. From this, we would state that 
the convergence of our method is excellent for stars with 
[N/Fe]> -0.5. 



4 RESULTS 

For each star the following parameters were determined: 
T eff , logy, [Fe/H], [C/Fe], [N/Fe] and [Ba/Fe] (Table g}. 
The CMD position of the stars is shown in Figure [l] There 
were 221 stars that were analysed with 14 having two spec- 
tra. An estimate of the uncertainty in the abundances de- 
termined was made using the doubly observed stars. The 
average standard deviation in the difference in [Fe/H] was 
±0.2 dex, [C/Fe] was ±0.2 dex, [N/Fe] was ±0.3 dex and 
[Ba/Fe] was ±0.6 dex. Among the 221 stars previously men- 
tioned there were nine stars in common between |M±2012"1 
and vL±07| (Table [5| . These stars were used for direct com- 
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Table 3. Results of the delta analysis performed on four stars to understand the effect of changing the 
input temperature by +100 K and the gravity by +0.3 dex. 



LEID 




[Fe/H] 




[C/Fe] 


[N/Fe] 






AT cS = +100 K A log g = 


+0.3 dex 


+100 K +0.3 dex 


+100 K +0.3 dex 


25062 


3965 


-0.1 


+0.1 


-0.3 -0.2 


+0.1 +0.3 


61085 


3965 


+0.2 


0.0 


-0.2 0.0 


-0.4 +0.1 


30013 


4485 


-0.1 


0.0 


0.0 0.0 


-0.1 +0.1 


57114 


4470 


+0.2 


-0.1 


0.0 0.0 


-0.1 +0.1 



Table 4. Stellar parameters determined for the 221 stars. The LEID is the star label given by |vL+07] Abundances from 
other studies are indicated. The offsets described in Section [4] and Figure [4] have been applied to [C/Fe] and [N/Fe]. The 
group number is described in Section [5] 



LEID T cff logs [Fe/H] [C/Fe] [N/Fe] [O/Fcf [Na/Fe] a [Ba/Fe] S(3839) 6 CH(4300) 6 Group 



15022 


4400 


0.9 


-1.8 


-0.5 


0.3 


0.44 


-0.36 


0.1 


0.108 


0.229 


2 


15023 


4450 


1.0 


-1.6 


-0.6 


0.1 


0.42 


-0.29 


-0.2 


0.075 


0.264 


1 


16019 


4800 


1.6 


-1.7 


0.7 


-0.1 


0.64 


0.01 


0.8 


0.475 


0.400 


2 


16027 


4750 


1.7 


-2.3 


-0.1 


-1.0 


0.46 


-0.16 


-0.1 


0.089 


0.239 


1 


17029 


4800 


1.6 


-1.9 


-0.8 


-0.9 


0.35 


0.49 


-0.2 


0.090 


0.135 


1 


17046 


4850 


1.7 


-2.3 


-0.6 


1.0 


0.14 


0.16 


1.0 


0.049 


0.120 


3 


37318 


3950 


0.8 


-1.0 


-0.4 


2.0 


0.35 


1.03 


0.5 


0.411 


0.337 


1 


38115 


4150 


0.8 


-1.6 


0.4 


0.1 


0.28 


0.14 


0.2 


0.397 


0.422 


2 


44449 


4100 


0.6 


-1.5 


-0.6 


1.5 


-0.73 


0.47 


0.5 


0.572 


0.409 


3 



" IJfcPIO 
HvL+07 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance 
regarding its form and content.) 



-1.2 
[Fe/H] (J&P10) 



Figure 3. Direct comparison of the metallicities found for each 
star in this study and the abundances found for the same stars 
by |JfcP10| (with the one-to-one line shown) 



parison of the [C/Fe] and [N/Fe] abundances as J&P10 did 
not determine these elements. 



The abundances determined were compared to |JfcP10| 
and |M+2012| directly and indirectly. Plotting the metallicity 
found by this work against that of |JfeP10| shows a good cor- 
relation (Figure [3}. Most of the stars had [Fe/H] from -2.0 
to —1.5, resulting in a large concentration around these val- 
ues. As would be expected when results derived from spectra 
with a resolution of 1,600 is compared to spectra with a res- 
olution of 20,000, there is a spread of values. About 70% of 
the stars were within ±0.2 dex of the value found by |JfcP10| 



The average difference for all the stars was +0.1 dex: |JfcP10| 
found a slightly more metal-rich value than that found by 
this work. 

The nine stars in common between this study and 
M+2012| (Table [5} had a range of T cff from 4050 to 4400 K 
and in metallicity from — 1 > [Fe/H] > —2. According to the 
abundances determined by |M+2012| four would be classed 
as CN-strong ([N/Fe]> 0.6), while four were CN-weak with 
higher carbon abundances. This meant that they provided 
a good sample across the range of stellar parameters which 
was not clumped at a particular metallicity or [X/Fe]. 

A comparison of our [C/Fe] and [N/Fe] to that of 
(20121 showed that there was a systematic 



Marino et al. 



offset between the two datasets. In order to place this work 
onto the same abundance scale as M+2012 it was necessary 
to apply some systematic corrections to our values. Using 
the eight of the nine stars in common (LEID 36179 was ex- 
cluded as an outlier), it was found that the [C/Fe] needed the 
addition of +0.45 dex, while the offset in [N/Fe] was —0.55 
dex. This offset was likely caused by continuum placement 
and the positive feedback between [C/Fe] and [N/Fe]. If the 
[C/Fe] is too low, then the [N/Fe] will have to be increased 
in order to reach the same strength of the CN region. The 
comparison of the corrected [C/Fe] and [N/Fe] values are 
shown in Figure [4] 

The significant outlier in Figure |4]is LEID36179. Fi gurc 
[5] shows the synthetic spectra produced using our parame- 
ters and the parameters of M+2012 for this star. This clearly 
shows that the |M+2012| C and N abundances do not provide 
a good fit to the CN and CH spectral features in the vL+07| 
spectrum. One note about this star is that [M+2012 deter- 
mined an [0/Fe]= -0.1 while |JfePl5"l have [0/Fe]= -0.42. 
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Table 5. Stellar parameters determined (this work = SCW) for the nine stars in common with |M+2012| The LEID is the 
star label given by |vL+07| and the ID is the numbering scheme from |M+20"T2| The [Fe/H], [C/Fe], [N/Fe] from |M+2012] 
arc shown for comparison. 



LEID 



ID 



vL+07 M+2012 



[Fe/Hl 



log 3 SCW M+2012 



[C/Fe 



SCW 



M+2012 



[N/Fe 



SCW 



[O/Fe 



M+2012 




J&P10 




M+2012 



35090 


246585 


4100 


0.8 


-1.3 


-1.64 


-0.55 


-0.65 


1.05 


1.55 


-0.32 


-0.25 


36156 


245724 


4150 


0.7 


-2.0 


-1.97 


-0.45 


-0.44 


-0.05 


0.20 


0.43 


0.30 


36179 


244812 


4050 


0.8 


-1.5 


-1.40 


-0.75 


-0.17 


2.05 


1.40 


-0.42 


-0.10 


37253 


242745 


4150 


0.7 


-1.9 


-1.91 


-0.25 


-0.35 


0.15 


0.10 


0.43 


0.37 


38115 


241359 


4200 


0.8 


-1.6 


-1.69 


0.35 


0.17 


0.05 


0.26 


0.28 


0.40 


46150 


224500 


3950 


0.6 


-1.6 


-1.51 


0.25 


0.40 


0.65 


0.20 


0.44 


0.49 


46194 


225246 


4200 


0.9 


-1.9 


-1.87 


-0.65 


-0.61 


0.65 


0.70 


-0.05 


0.19 


48235 


220325 


4050 


0.7 


-1.5 


-1.61 


-0.65 


-0.45 


1.15 


1.20 


-0.32 


0.00 


51091 


215367 


4400 


1.1 


-2.0 


-1.93 


-0.55 


-0.98 


0.75 


0.90 


0.28 


0.15 



i 

I 1.2 



ta) 








® 




















♦ 






♦ 






♦ 




















♦ 










♦ 











0.6 1.0 
[N/Fe] (M+2012) 




-0.4 -0.0 
[C/Fe] (M+2012) 

Figure 4. (a) [N/Fc] and (b) [C/Fe] abundances found by this 
work and |M+2012] The dashed line is the one-to-one line. An 
offset of —0.55 dex has been applied to our initial [N/Fe] abun- 
dances and +0.45 dex applied to our initial [C/Fe] abundances. 
This was done to place our results on the same abundance system 
as [M+2012] The circled star is LEID36179, the outlier discussed 
in Section |4] 




0.2 



LEID36179 



M + 2012 
Synthetic 
Observed 



4100 4150 4200 o 4250 
Wavelength (A) 



4300 



Figure 5. LEID36179 is the significant outlier on figure |4p,. We 
determined [N/Fe]=2.35, while [M+2012] determined 1.4 dex. The 
thick solid line shows the synthetic spectrum that results using 
|M+2012| s abundances adjusted using the offsets described in Sec- 
tion |4 and Figure [4] The thin solid line is the observed spectrum 



n v 



from |vL+07| and the dotted line is our fit to this spectrum. 



This is the largest difference for these nine stars in terms of 
the oxygen abundance determined by the two groups. The 
overall results in this research are not compromised by this 
outlier. 

A qualitative comparison (Figure [6} of the [C/Fe] to 
[N/Fe] shows that the automated spectral matching method 
worked well: it did not find stars that were both carbon and 
nitrogen rich, the upper-right quadrant on that plot. There 
was not the clear bifurcation that was seen by |M+2012l but 
there was clearly the same anti-correlation between [C/Fe] 
and [N/Fe]. Of interest were the very nitrogen-rich, carbon- 
poor stars that are outside the range found by |M+2012| 
These could be related to the way the spectral matching 
worked. It kept searching for the parameters that minimized 
the x 2 j which could sometimes lead to extremely large or 
small values as the \ 2 difference between the steps in [C/Fe] 
or [N/Fe] were very small. This was the cause of the ex- 
tremely low [N /Fe] values that are not shown on Figure [6] 
These stars have so little detectable nitrogen that that the 
automated process determined values as low as [N/Fc]= —3. 
This combination of the parameters and the resolution gives 



© 2012 RAS, MNRAS 000,[T}{T7] 



8 J. D. Simpson, P. L. Cottrell and C. C. Worley 



Si.o 



+ #■♦ + ♦. ♦ 
+ +♦ + * 

♦ ♦ "t ♦ 

♦hu- ♦ ♦ ♦ ♦ ♦ 

♦ ♦ +• ♦ 



+ + + « 



+ « . .4- . ♦+♦+*♦♦ . + 

+ :+-!-: : * : : <• +:+• + 



[C/Fe] 



Figure 6. The [C /Fe] and [N/Fe] for this work (filled diamond) 
and the |M+2012] (+). The abundances were only determined to a 
precision of 0.1 dex due to the spectral resolution of the data. It 
shows the expected anticorrelation, as observed by |M+2012| This 
plot excludes ten stars that had extremely low [N/Fe] abundances. 



just an upper limit on the [N/Fe], rather than a definitive 
abundance. 



5 CLUSTERING ANALYSIS 

fc-means clustering is a method to create k groups from a 
dataset of n entries using m parameters. It minimizes the 
distance of each member from the group (or cluster) centre. 
In the work presented here, the method was implemented 



using the R statistical package ( R Development Core Team 
|2011[ ), a freely available system for statistical computation 
and graphics. The default R algorithm of jHartigan fc Wong| 
( 1979 1 was employed in this work. 

There are two caveats for fe-means clustering. First, 
that the Euclidean distance is used for determining the clus- 
tering. This could require the observational data be trans- 
formed to give an appropriate range of values. This was not 
done in this work as the range of values was small and of 
similar magnitudes. Second, the number of groupings must 
be appropriate for the dataset. Too few or too many group- 
ings can lead to less definitive results. In this case, there was 
experimentation with a range of the number of clusters. 

|Gratton et al" (20111 used fc-means clustering on the 
| JfcP10| dataset to create seven groupings of stars in u> Cen. 
They used the [Fe/H], [a/Fe], [Na/O] and [La/Fe] to show 
there are three metallicity regimes in the cluster with a 
range of a, light element and neutron capture properties. 
Their analysis confirmed each metallicity regime has its own 
correlations and anticorrelations amongst the different ele- 
ments. Of particular interest was the very metal-rich group 
that showed a Na-0 correlation and had the highest aver- 
age neutron-capture abundances of the groups in the cluster 
they identified. 

We recognized that our low resolution data had larger 
uncertainties in each parameter than J&P10. This would 
make group identification harder. In order to first under- 
stand the groups that could be found when carbon and nitro- 



gen were used in conjunction with some of the abundances 
that Gratton et al.|(|2011|) u sed, we decided to combine the 
Marino et al.| ( |2011[ |2012| ) datasets. There were 74 stars 
which had an [Fe/H] and [C,N,0,Na,Ba/Fe]. These six pa- 
rameters were chosen to look for groupings that could be 
also found in our dataset of 221 stars. Neither Marino et al. 



(20111 nor M+2012 determined an a-element abundance. 



However, Gratton et al. (2011) found that there was a posi- 



tive correlation between metallicity and their a-element av- 
erage for their groups. 

In our analysis of |M+2012| it was found that a clus- 
tering analysis of five groups gave the best results. Two 
groups would simply split between Na-strong and Na-weak 
stars. Four groups did not isolate the stars with extreme 
sodium abundances, with these stars being found in a group 
that consisted of a large range of metallicities and oxygen 
abundances but the same barium abundances. Although five 
groups were used, two of the groups were very similar in their 
properties so were combined into one group (#3a & #3b). 

The clustering analysis was then performed on our re- 



sults using [Na,0/Fe] from |JfeP10| and [C,N,Ba/Fe] and 
[Fe/H] from this work. Although we were able to define five 
groups from the Marino et al. (2011 2012| ) data, only four 
groups were attempted with our data. This was because one 



of those groups (group #3b) found in Marino et al. (2011 



2012 1 was primarily the result of it having higher barium 
abundance for its metallicity. This level of precision was 
not available in our data. As such, the same four overall 
groups were found with group #3 having no subset that 



was s-process enhanced. Although Gratton et al. (20111 de- 
termined seven groups from the |JfcP10| dataset, we did not 
think that our smaller sample with lower precision abun- 
dances would allow us to define that many groups. Any fur- 
ther division created signal out of noise. 

The cluster centres of each group are given in Tables [6] 
and[7| In the context of this paper, group #1 will be referred 
to as metal-poor, groups #2 & 3 as intermediate metallicity 
and #4 as metal-rich. By way of comparison, our group #1 
is the equivalent of groups 4 & 6 in |Gratton et al.| ( |2011[ ) , 
our #2 is their 1 & 5, our #3 is their 2a & 3, and our #4 is 
their 2b. 



5.1 Groupings on the colour-magnitude diagram 



5.1.1 Marino et al. , '2011 2012) 



The stars of 



Bellini et al 



Marino et al. (2011 2012 I were combined with 



( 2009 1 and its high precision photometry. Fig- 
ure |7| shows "the V-(U-V) and U-(U-V) CMD with the 
five groups plotted. U photometry is used here instead of B 
as it was found to create a greater separation between the 
groups. Not all of the stars in our 221 star sample have U 
photometry but there was no bias in the stars which did, 
and did not, have U magnitudes determined by |Bellini et al.| 



(2009) 



The metal-poor group is clearly on the hot-edge of the 
giant branch with little colour spread. As the metallicity 
increases, the stars are found at cooler positions. There are 
four stars in group #4 which are likely members of the RGB- 
a ( Pancino et al.|2000 l, the most metal-rich sequence in the 
cluster. These stars are found to be the most enriched in all 
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Table 6. The cluster centres found using fc-means clustering analysis of the 221 stars from the |vL+07| spectral library. 
For each parameter, the average and standard deviation are given, along with the number of stars in each group. The 
symbol is the identifier used for that group in the figures. 



Group 


Symbol 


N. stars 


[Fe/H] 


[C/Fe] 


[N/Fe] 


[O/Fc] 


[Na/Fe] 


[Ba/Fc] 


1 


* 


62 


-1.8±0.3 


-0.4 ±0.2 


0.3 ±0.3 


0.3 ±0.1 


-0.1 ± 0.3 


-0.4 ±0.5 


2 


□ 


86 


-1.7 ± 0.2 


0.0 ±0.3 


0.4 ±0.3 


0.4 ±0.1 


0.0 ±0.2 


0.4 ±0.3 


3 


▲ 


49 


-1.7 ±0.3 


-0.5 ±0.2 


1.0 ±0.3 


-0.3 ±0.4 


0.3 ±0.2 


0.5 ± 0.4 


■1 


X 


24 


-1.3 ±0.4 


-0.6 ±0.2 


1.8 ±0.4 


-0.2 ± 0.2 


0.7 ± 0.2 


0.8 ±0.4 



Table 7. The cluster centres for the 74 stars from the |Marino et al| ( |2011||2012| l data. The column descriptions are the 
same as those in Table [6] 



Group Symbol N. stars [Fe/H] [C/Fe] [N/Fe] [O/Fe] [Na/Fe] [Ba/Fc] 



1 


* 


17 


-1.9±0.1 


-0.5 ±0.3 


0.5 ±0.3 


0.3 ±0.1 


0.1 ± 0.2 


-0.4 ± 0.2 


2 


□ 


18 


-1.6±0.2 


0.3 ±0.2 


0.4 ± 0.2 


0.6 ±0.2 


0.1 ± 0.1 


0.4 ± 0.2 


3 


A 


19 


-1.6±0.1 


-0.6 ±0.2 


1.3 ±0.2 


0.0 ±0.2 


0.4 ±0.2 


0.3 ±0.2 


3b 


A 


13 


-1.5 ±0.1 


-0.5 ±0.2 


1.4 ±0.1 


-0.1 ± 0.2 


0.6 ±0.1 


0.9 ±0.2 


4 


X 


7 


-1.1 ± 0.1 


-0.3 ±0.3 


1.6 ±0.1 


0.2 ± 0.2 


1.0 ±0.1 


0.8 ±0.1 




Figure 7. (a) V-(U — V) and (b) U-(U — V) colour-magnitude 
diagrams of the M+2012 stars, with different symbols for the 
different groups as described in Table [7] The small dots are the 
whole Bellini ct al. ( 2009]l dataset. The four stars with the faintest 
U magnitude are part of the metal-rich group and are on the 
RGB-a. The AGB star described in Section l5.1.1l is circled. 



the elements measured, except sodium where all of group 
#4 are equally enhanced. 

As would be expected from their metallicity ranges, 
groups #2 and #3 have very similar magnitude and colour 
ranges, while group #3b is slightly cooler with a tight locus 
like group #1. As discussed in the following sections, group 
#3b could be an evolutionary extreme of group #3 or the 
end point of the sodium-oxygen anticorrelation trend before 
the stars begin to be sodium-oxygen correlated. 

The correlation between metallicity and colour is ex- 
pected from stellar evolutionary models. The small spread 
in [Fe/H] in group #1 correlates with a small spread on 
the CMD. Conversely, the larger spreads in the interme- 
diate metallicity groups are shown on the CMD. On the 
Tefi-log g plane, a straight line was fitted through the group 
#1 stars, which could be used to produce a vertical plot of 
temperature deviation from the expected temperature for 
that gravity. The standard deviation of the temperatures 
in group #1 using this technique was 30K. Group #2 were 
100 ± 75A" cooler, #3 were 40 ± 90A", #3b were 80 ± 30K 
and #4 were 300 ± 60K. 



It is possible that star ID 215367 (Marino et al. 2011 



numbering as this star is not in the |vL±07| library) is an 
AGB star (circled in Figure f7|, as it sits slightly to the left 
of the main locus of group #1 stars on Figure[7] This is most 
obvious in a (U,U — V) CMD (Figure (7f}), as this region of 
the GB is at an almost constant U magnitude. Here this star 
is at a much brighter U magnitude than would be expected 
for its U — V. This star is found to have the lowest [C/Fe] 
(and [C/H]) of all the stars measured (discussed in Section 
5.3 and shown in Figure 10 1. It is the most barium-rich of 



all the group #1 stars (Section 5.4 1 



5.1.2 This work 



Plotting the colour-magnitude diagra ms (Figure p]) f or our 



results has the same findings as with |Marino et al. (2011 



2012|. Compared to Figure [7] there is more blending of the 



groups which is caused by the larger standard deviations in 
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5.2 Sodium-oxygen ant icor relation 



Figure 8. The same as Figure[7]but using the 221 stars from this 
study. Symbols are the same as in Table [6] The sequences are not 
quite as neatly defined as with the |Marino et al~| (2011 2012J data 
but do show the same effects. Wc find that the metal-poor group 
(#1) is concentrated on the hot side of the giant branch and the 
metal-rich stars (group #4) are found in the anomalous RGB 
sequences. There are a number of potential and probable AGB 
stars found. The nine stars that are between 13.0 > V > 13.5 
and (U—V) < 1.4 are highly likely to be AGB stars (surrounded 
by a rectangle). Only 172 of 221 stars are plotted on this diagram 
as not all the stars had U photometry in the |Bellini et al.| ( |2009| ) 
photometric library. 



the average metallicities of the groups than for |Marino et al.| 
(2011 2012 I. However at the bright end of the giant branch, 



there is still a tight locus of the most metal-poor stars. The 
metal-rich stars are also the coolest as would be expected 
and the two intermediate metallicity groups are inseparable. 

Our larger sample size has found nine stars that are in 
a well-defined AGB sequence. These stars come from the 
metal-poor and intermediate metallicity groups and are all 
found in the carbon- and nitrogen-poor quadrant of the CN 
diagram. These and other AGB stars will discussed in Sec- 
tion [6] in the context of the potential absence of CN-strong 
stars on the AGB. The J&P10 dataset did not extend to 
magnitudes fainter than V = 13.5, which means there are 
only two members of the RGB-a metal-rich giant branch 
( Pancino et al.|2000 l. 



The fc-means clustering of Marino et al.| ( |2011| 2012 1 dataset 
split the Na-0 anticorrelation at about [Na/Fe]= 0.3, with 
groups #1 & #2 being oxygen-rich, sodium-poor and groups 
#3 & #4 having the opposite properties (Figure [9} . The 
results broadly confirm the results found by |JfcP10| and 
|M+2012| They both found that the intermediate metallic- 
ity stars displayed the well-known anticorrelation of mono- 
metallic clusters, while the metal-poor and -rich groups do 
not share this bimodal behaviour. The same result is ob- 
served for the clustering analysis of our data. 

Plotting the abundances with respect to hydrogen in- 
stead of iron (Figure [9J;), the group #1 stars have the same 
[O/H] as group #3 with [Na/H] different by ~ 1 dex. Group 
#2 stars are about 0.5 dex enhanced in [O/H] compared to 
group #1 stars and potentially slightly enhanced in sodium. 
If group #2 formed from the the pollution of group #1 stars, 
this gas and dust would be enhanced in overall metallicity, 
oxygen and sodium. However if group #3 were formed after 
group #1, this would require an increase in metallicity and 
sodium and a depletion in oxygen. 



D'Antona et al. (2011 I would have groups #3 and #3b 



at opposite ends of the evolutionary sequence that forms the 
sodium-oxygen anticorrelation. They theorized that the first 
stars of the second generation were sodium-rich, oxygen- 
poor and formed from pure AGB ejecta. As this ejected ma- 
terial was diluted by mixing with the cluster's surrounding 
interstellar medium, progressively oxygen-richer, sodium- 
poorer stars would form. This would continue until the AGB 
ejecta from these stars becomes the majority contributor to 
the gas of the cluster and produces sodium-rich stars. In this 
way, #3b could represent an extreme limit to sodium en- 
hancement. The extension of group #2 into extremely high 



oxygen values is seen in the Marino et al. (20111 dataset, 
but is not evident to the same extent in the larger |JfcP10| 
dataset. 



Both Marino et al. (2011) and our study show the cor- 



relation of sodium and oxygen in the most sodium-rich 
stars. This feature is not seen in mono-metallic clusters. 
The metal-rich stars of group #4 form this sequence, which 
would imply that they are the youngest stars of this sam- 
ple, forming from the ejecta of the last supernovae in the 
cluster. Gratton et al. ( |2011[ ) found this group to be the 
most enhanced in a-elements as well. The correlation is most 
dramatically seen in Figure[9^ which used the J&P10 abun- 
dances. Here the correlation extends from [0/Fe]= —0.6 to 
+0.6. 



5.3 Carbon and nitrogen 

As with sodium-oxygen, there are roughly three regions of 



chemical abundance for carbon and nitrogen (Figure 10 1: #1 



are C- and N-weak; #2 are C-strong, N-weak; and #3 and 
#4 are C-weak, N-strong. There is a range of abundances in 
each grouping with #1 showing a range of 1 dex in [N/Fe] 
and about 0.6 dex in [C/Fe]. The intermediate metallicity 
groups have a range of 0.6 dex in [N/Fe] and [C/Fe]. Group 
#4, which has very few stars, ranges from [C/Fe]= —0.65 to 
+0.3. 

Groups #3 and #3b have similar ranges of [C/Fe] (Fig- 
ure 10 a) but the #3b group form an upper limit for ni- 
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Figure 9. The sodium-oxygen (anti)correlation for u Cen of (a) 
this work and (b, c) |Marino et al.| ( pOll] [2012 1. Symbols as de- 
scribed in Tablcs[6]&|7|and uncertainties are shown in the bottom- 
left of each diagram. The small dots in (b) & (c) are the fuller 
|Marino et al.| ( |201 1^ dataset. They illustrate that features ob- 
served (for instance the correlations within group #2 and #4) 
are not due to small numbers of stars. The sodium-oxygen anti- 
correlation is clearly present between groups #2 and #3. There 
is little separation between the two subgroups of #3, though po- 
tentially #3b is slightly more sodium-rich than #3. Group #1 is 
found to be the most sodium-poor on average but with an inter- 
mediate oxygen, both with respect to iron and hydrogen. 
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Figure 10. Nitrogen against carbon for (a) this work and (b, c) 
|Marino et aTT1 < |201 1|[2012] ) . Uncertainties are shown in the bottom- 
left of each diagram. The C-N plane looks broadly similar to the 
Na-O plane. The metal-poor group (#1) is found to be deficient 
in both carbon and nitrogen, while the intermediate metallicity 
groups are bimodal. Group #3b shows an extreme enhancement 
of nitrogen. When considered with respect to iron (a, b), group 
#4 stars are very similar to those of group #3 in terms of their 
carbon abundance. However when they are considered with re- 
spect to hydrogen (c), their high metallicities result in a separate 
sequence, with some up to 0.8 dex higher than the solar value 
in nitrogen. The potential AGB star (Se ction |5.1.1| is the most 
carbon-deficient star in the |Marino et aiT| | |2011||2012[ ) sample (b, 
c). 
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Figure 11. The total C+N+O with respect to iron for the 
ct al. (2011 2012) dataset (using solar abundances from 
& Grcvcssc 1989). Symbols as described in TableJT] Uncertainties 
are shown in the bottom-right of each diagram, (a) the oxygen- 
rich groups (#1 & #2) have oxygen as their dominant CNO 
species, (b) while in the nitrogen-rich groups (#3 & #4) it is 
nitrogen. Unlike mono-metallic clusters, the total C+N+O is not 
constant. 



trogen. This is reinforced when the stars are plotted with 



[C/H] versus [N/H] (Figure [10};) . The group #3b stars have 
an average [N/H]= 0.01 ± 0.06 suggesting a maximum pos- 
sible pollution for stars of this generation. Group #4 stars 
are found to have even higher nitrogen abundances with the 
four RGB-a ( |Pancino et al.|2000[ ) stars having [N/H]> 0.45 
([N/Fe]> 1.4). 

In mono-metallic clusters, the total C+N+O has been 



found to be constant (e.g. Carretta et al. 2005 1. In the case of 



to Cen there is a large spread (Figure 11 1. It was found that 
groups #1 and #2 have a straight-line relationship between 
[O/Fe] and [C+N+O/Fe] but with the other two groups hav- 
ing no correlation. Between [N/Fe] and [C+N+O/Fe], it is 
instead groups #3 and #4 that have the positive correlation. 
This illustrates that in the first two groups, it is nitrogen 
which is the dominant CNO species, while in the other two 
groups it is oxygen. This is what would be expected based 
upon their relative abundances of nitrogen and oxygen. 



In |Marino et al.| ( |2011[|2012[ ) (see our Table[7| there are 
39 (53%) stars in the N-rich groups and 35 (48%) in the N- 
poor group, while in our work (see our Table[6]) there are 148 
(67%) N-rich stars and 73 (33%) N-poor stars. We attribute 
this to a selection effect. The split in nitrogen groupings 



maps to the split in oxygen. In the Gratton et al. (2011) 



analysis, they had 574 (72%) stars in O-poor groups (N- 
rich) and 223 (28%) stars in O-rich groups (N-poor). So our 
groupings follow that of Gratton et al. (2011) in terms of 



the split of N-rich and N-poor stars. Our selection of stars 
was effectively based upon the selection of |JfcP10] through 
the use of their [Ca/Fe] and [O/Fe] abundances. 



5.4 



s-process 



A large range of s-process abundances is not typically ob- 
served in mono-metallic clusters. However, ui Cen does ex- 
hibit a very large range of abundances in elements such 
as barium. With a simple correlation between [Fe/H] and 
[Ba/Fe], the groupings were as would be expected from pre- 
vious work. The metal-poor stars all had sub-solar abun- 
dances of barium with respect to iron, but with a very 
large range. Some stars had [Ba/Fe]= —0.8 while others had 
greater than the solar abundance in [Ba/H]. 

Unlike the other element combinations investigated in 
this study (C-N & Na-O), barium does not exhibit any dif- 
ferentiation between the intermediate metallicity groups of 
#2 and #3. This could require that any enhancements in 
metallicity and s-process happened before the anticorrela- 
tions. 

Here the separation of group #3 from #3b was most 
evident (Figure [T2[ a). The two groups had about the same 
metallicity range but a difference of about 0.5 dex in [Ba/Fe] . 



A caveat to this is that in Marino et al. ( 2011 1 there were few 



stars at this metallicity/s-process abundance combination. 
They comment that at high barium abundances the lines 
get too strong for accurate measurement. However we have 
kept them as a separate grouping due to their separation in 
the other abundance planes that we investigated. 

Neither |Marino et al.| (|2011| [20121 nor this work de- 



termined any r-process abundances. |JfcP10| did determine 
[Eu/Fe]. As noted by Gratton et al. (2011 1 there is no trend 



of [Eu/Fe] with [La/Fe] in J&P10 indicating that the contri- 
bution of the r-process to neutron-capture elemental abun- 
dances in the cluster are negligible. It would appear that 
there is a potential increase in [Eu/H] at high [La/H] and 
that the [La/Eu] abundance is dominated by changes in 
La. There is a larger scatter in [La/Fe] in the the metal- 
intermediate and -rich groups (almost 1 dex) but overall, in 
these groups, this ratio is constant. 



6 LACK OF CN-STRONG AGB STARS 



Campbell et al. ( 2010 1 found evidence that in some globular 



clusters, there are no CN-strong stars on the AGB. This was 
first found by |Norris et all ( [T98l| in NGC 6752. This lack 
of CN-strong stars on the AGB is in contrast to the RGB 
where the stars are equally CN-weak and CN-strong in all 
clusters, including lj Cen. In the case of |M+2012| there is 
only one potential AGB star. This star has a low [N/Fe] 
and [C/Fe], which would fit with the findings of Campbell] 
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Figure 12. The correlation between the metallicity of the stars 
and their s-process abundance for (a) this work and (b, c) |Marino| 
|et al.| | [2uTT||2012^ . Uncertainties are shown in the bottom-right 
of each diagram. Only group #4 docs not form part of the cor- 
relation, with the s-process abundance constant with respect to 
iron. One of the causes of the separation of group #3b from group 
#3 is the roughly 0.5 dex difference in their [Ba/Fc]. 
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Figure 13. (a) The four groups with th eir CN [S(3839)] and CH 
[CH(4300)] indices as found by |vL+07| T he large diamonds are 
the AGB stars identified in Section |5.1.2| The curved line sepa- 
rates the oxygen-rich groups, defined as [0/Fe]> 0.17 (below the 
line) from the oxygen-rich stars. This serves as a divider between 
CN-weak and CN-strong stars. The vertical dashed line divides 
the HB from the GB in these diagrams. The small dots in all sub- 
plots show all the stars in the |vL -)-07| dataset. (b) Highlighting 
all the 55 AGB stars identified in t he entire |vL+07| libr ary based 
upon their position in the CMD of |Bellini et al.| p009\ . 



|et al.| ( |2010p . From Section |5.1.2| th ere are nine potential 
AGB stars in the subset of vL+07 for which [C/Fe] and 
[N/Fe] were determined. 

|vL+07| measured two indices of interest: the S(3839) 
and CH(4300) indices^ The first measured the CN band 
around 3839A and the secon d CH band at 4300A (the G 
band). Using the definition of Harbeck et al. (20031, where 
F\i-\2 is the average flux level between the two wavelengths 
in angstroms: 



S(3839) = -2.5 log 
and 

CH(4300) = -2.5 log 



F3861-3884 
F3894-3910 



2F4285-4315 
F4240-4280 + F4390-446O 



(1) 



(2) 



Figure 13 1 shows the S(3839) and CH(4300) indices 



3 It should be noted that these values used |vL+07| normalized 
spectra. 
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strength for all the stars determined, with a curved line di- 
viding the O-rich stars ([0/Fe]> 0.17) (Figure [9| from the 
O-poor stars. This line takes the form of a parabola. The 
clump of stars in the lower-left quadrant of both subfigures 



of Figure 13 with CH(4300)< 0.9 are HB stars. This has 



been separated from the GB by a vertical dashed line. The 
groups of stars which were deficient in nitrogen (see Figure 
10 1 are almost all below the dividing line, while group #3 (C- 
weak, N-strong) stars are placed above the line. The (metal- 
rich) group #4 are found to be on both sides of the line. 
For this group there is a much larger variation in S(3839) 
than in CH(4300), which would correspond to their larger 
range of [N/Fe] than [C/Fe]. The nine AGB stars have small 
to medium strength CN and CH bands based upon their 
strengths of their indices. 

All this gives good evidence that we can use the S(3839) 
and CH(4300) indices of |vL+07| to define the CN strength 
of the stars. To extend our sample of AGB stars from just 
those for which a [C/Fe] and [N/Fe] was determined, the 
photometry of |Bellini et al.|p0 09) was used to separate the 
AGB from the RGB from 12.6 <V< 14.2. By eye, over 200 



AGB stars were defined from the Bellini et al. ( 2009 1 sample 



(out of over 60,000 probable cluster members). These stars 
then gave 55 AGB stars in the overall |vL+07| sample for 
which there was also Bellini et al. ( 2009 1 photometry. 



All the |vL+07| stars identified as AGB stars have low 
S(3839) and CH(4300) indices (Figure pja). In the case of 
the CN index, they are all less than 30% of the maximum 
value measured for any star by |vL+07| The vertical line 
places 18 of the 55 stars (33%) in the HB part of the diagram. 
Some of these stars may not be HB stars but are too hot 
to easily determine their CN strength from these indices. 
Of the remaining 37, only two of these AGB stars (5%) are 
significantly above the line. From this we would infer that 
more than 90% of this AGB are CN-weak. 

It should be noted that in the case of uj Cen, this refers 
to one particular AGB. Due to the range of ages, metallic- 
ities and helium content, there will be potentially as many 
distinct AGBs as there are distinct RGBs. This means that 
it is almost certain that some of the stars that are being 
referred to here, and in other works as RGB stars, will in 
fact be AGB stars. This level of contamination will also be 
affected by how many of stars actually ascend the AGB. 



D'Cruz et al. (20001 estimated that as many as 30% of the 
horizontal branch stars are extreme in their colour, which 
could mean that they evolve directly to white dwarfs, never 
passing along the AGB. 

In the case of the AGB identified in this work, the as- 
sumption would be that it is the AGB is associated with the 
hottest and most metal-poor RGB. This RGB is thought 
to be not helium-rich, being associated with the red main 
sequence. A low helium abundance would not support the 



theory of Norris et al. ( 1981 1 that the lack of CN-strong 



stars was the result of them being helium-rich and therefore 
evolving directly from the HB to the white dwarf phase. 



7 DISCUSSION 

The clustering analysis identified four groups of stars, which 
would have had multiple star formation events within them. 



Each grouping showed a range of abundances in every ele- 
ment that was incorporated into the fc-means clustering. 

Our primordial population is group #1. These stars are 
deficient in every element except oxygen. This implies that 
the material from which the cluster formed had not been 
polluted with material from AGB stars, if it is assumed that 
all generations of the cluster formed in the same spatial lo- 
cation and not as the result of a merger. Within this first 
generation there is still a metallicity range, greater by an 
order of magnitude than what is observed in mono-metallic 
globular clusters. 



D'Antona et al. (2011 1 used the yields of massive AGB 



stars to model the sodium-oxygen anticorrelations (and cor- 
relation in the case of the very metal-rich stars) in uj Cen. 
Using the 300 stars of Marino et al. (20111 they identified 
that for all the stars in the metallicity range of —1.9 < 
[Fe/H] < —1.3, there was a Na-0 anticorrelation, as ob- 
served in our groups #2 & #3. This implies that the process 
through which the anticorrelations formed will be broadly 
similar to that which formed the anticorrelations in mono- 
metallic clusters. 

The most metal-rich stars in the cluster ([Fe/H]> —1.3) 
show a direct sodium-oxygen correlation (Figure Rjl. This, 



D'Antona et al. (20111 suggested, would be the result of 



stars forming from the pure ejecta of super- AGBs. In these 
stars the sodium and oxygen yields are correlated with the 
yields of both elements decreasing with stellar mass to about 
6.5 Mq, then the yields begin to increase again. This is 
caused by the changing temperature of the convective enve- 
lope in these stars. By the time the most metal-rich stars 
in the cluster form, this AGB material would not be di- 
luted as any residual ISM would have been blown away 
and/or consumed in preceding star formation episodes. This 
requires that future star-forming events would form from 
just the ejecta of AGB stars. Our analysis finds that these 
most metal-rich stars, with their high oxygen and sodium, 
are nitrogen-rich, but with no large enhancement of the s- 
process abundances found for their metal-poorer predeces- 
sors at high sodium but low oxygen. 

In figure 3 of D'Antona et al. (20111, there is a 



possible scenario for the production of sodium-oxygen 
(anti)correlation(s) in a mono-metallic cluster (in their case 
NGC 2808). Using the model of |D'Ercole elaT] ( [20081 [2010 \ 
they used super- AGB star yields to show a qualitative match 
to the Na-O anticorrelation and correlation of uj Cen. In this 
model, the first stars of the second generation are the most 
oxygen-poor, forming from pure AGB ejecta which has ac- 
cumulated at the centre of the cluster. These stars will be 
very helium-rich (corresponding to the blue main sequence). 
This pure ejecta is then diluted by "pristine" gas from the 
surrounding interstellar medium, which reduces its sodium 
abundance and increases its oxygen abundance. 

Of course such a model cannot be directly applied to 
uj Cen with its large metallicity range. The "pristine" gas 
that dilutes the pure AGB gas must be already enhanced 
in metallicity and s-process elements. The carbon-nitrogen 
anticorrelation also requires that this "pure" AGB gas be 
enhanced in nitrogen with no corresponding enhancement 
of carbon. Then the dilution with pristine gas will result in 
depletion in nitrogen and the enhancement in carbon. 

One of the most curious aspects of uj Cen is the cor- 
relation between metallicity and s-process abundance. Here 
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we found that both aspects of the Na-0 and C-N anticor- 
relation in the metal-intermediate groups showed the same 
range of s-process abundance. This would imply that any 
s-process enhancement that these intermediate metallicity, 
second generation stars experienced must have taken place 
before the second generation of star formation took place. 

This s-process correlation has always been difficult to 
reconcile with other aspects of the chemical evolution of the 
cluster. The generally accepted location of the formation 
of s-process elements is low-mass AGB stars (M< 3Mq). 
These stars require a timescale of up to a gigayear to 
evolve, while the source of the sodium-oxygen anticorre- 
lation, higher-mass AGBs (M> 5Mq), will have already 
evolved. As such there must be a substantial delay in the 
star formation that took place to create the second genera- 
tion of stars. And then a subsequent delay to form the most 
metal-rich stars, which again are even more enhanced in s- 
process. The most metal-rich stars are however not enhanced 
to the same extent, with [Ba/Fe] at the same values as the 
most-enhanced stars of the second generation. 

A clue to the formation of this cluster is that the 
oxygen-poor stars are centrally concentrated. [jfcPIO found 
that 91% of the stars with [O/Fe]^ were within 10 arcmin- 
utes of the cluster centre in their sample. This would be our 
groups #3 and #4. They inferred that the oxygen-poor stars 
were also from the blue main sequence (bMS), which is also 
the helium-rich sequence. So this implies that the oxygen- 
rich, helium-rich stars are also nitrogen-rich, carbon poor. 



This matches with the theory of D'Ercole et al. (2008 



20101 that the oxygen-poor, helium-poor stars were the first 
of the second generation stars to form and did so with a 
central concentration as this was where the cooling flows 
brought the pure AGB gas. Their theory would then require 
that the oxygen-rich stars that form later will have no radial 
gradient, instead forming throughout the cluster. 

An alternative theory was proposed in the recent paper 



of Herwig et al. (20121. They identified galactic plane pas- 



sages as mechanism for stripping the cluster of gas. Their 
evolutionary sequence would have our group #1 as the pri- 
mordial population, followed by #3 forming and then #4. 
The last population formed would be #2 from the gas of 
#1 after all any residual gas had been stripped by a galactic 
plane passage. 

One key element that requires measurement for a large 
number of stars in ui Cen is helium. Do our group #2 dif- 
fer greatly in their helium abundance from group #3 and 
do the group #3 subsets have differences in helium as well? 
The most recent determination of helium abundances was 



by King et al. (2012 I who used isochrones to determine that 



the bMS had a Y = 0.39 + 0.02. |Dupree et al.| ( |2011[ ) directly 
measured helium in giant branch stars and found that those 
stars with high sodium and aluminium tended to show he- 
lium, while those with [Na/Fe]< 0.2 did not show helium in 
their spectra. In terms of the data presented here, this would 
be groups #3 and #4 (and the most sodium-rich stars of 
groups #1 and #2) which would potentially show helium 
in their spectra. So the polluters that formed that gas and 
dust for group #3 would also be enhancing the helium abun- 
dance of these stars. This fits into the model hypothesized 



by D'Ercole et al. (2008 20101 



tra were then used to create synthetic CMD for a variety of 
bandpasses, one of which was the U and V used in this work. 
They noted that CNO variations will cause the greatest vari- 
ation in the U and B bands due to the molecular lines in this 
region of the spectrum. Their base isochrone had a metallic- 
ity of -1.62 (Z = 0.001, Y = 0.246) with an age of 12 Gyr. 
Three other abundance combinations were investigated: one 
enhanced in nitrogen and sodium by 1.8 and 0.8 dex respec- 
tively and depleted in carbon and oxygen by 0.6 and 0.8 dex 
(CNONal); one with the same enhancement in C, Na and 
O but only enhanced by 1.4 dex in nitrogen (CNONa2); and 
the same as the first but with Y = +0.4 (CNONal+He). 
The C+N+O mass fraction was twice the reference value in 
CNONal, while CNONa2 had the same value as the refer- 
ence. This is similar to the hypothesized situation in uj Cen: 
our group #3 has a larger C+N+O total than group #2 



(Figure 111 and from D'Ercole et al. (2008 20101 group 



#3 would be enhanced in helium. They found that in the 
U - (U-B) CMD, the CNONal+He and CNONa2 were 
indistinguishable for much of the RGB. This would concur 
with what we have found in this work: that groups #2 and 
#3 are blended on the CMD. But of course, their work used 
a mono-metallic isochrone and could not be fitted to the GB 
of w Cen. 



8 CONCLUSIONS 



The |vL+07| spectral library provides a large sample of spec- 
tra of giant branch and horizontal branch stars in the glob- 
ular cluster uj Cen. Combining this spectral data with high 
precision photometric data and previously determined abun- 
dances of these stars, we determined the [C/Fe] and [N/Fe] 
for over 200 new stars in the cluster. These results confirmed 
the results of |M+2012| that the low metallicity stars have 
both low carbon and nitrogen, while the intermediate groups 
show an anticorrelation. 

Combining this carbon and nitrogen data with sodium, 
oxygen and barium abundances, and using fc-means cluster- 
ing, we found four (potentially five) groupings of stars (with 
metallicity averages from Table [7|: 

• A low-metallicity group (([Fe/H]) = —1.9) which is defi- 
cient in carbon and nitrogen relative to iron, but not oxygen. 

• Two intermediate metallicity groups (([Fe/H]) = —1.6) 
which are anti-correlated in Na-O and C-N. Group #2 are 
high in C and O, while group #3 are high in Na and N. 
Both these groups have the same metallicity and barium 
range. There is a potential subgroup of group #3, which 
was labelled #3b which has a higher barium than the rest 
of the group and higher Na and N without any differences 
in its O and C. 

• A metal-rich group (([Fe/H]) = — 1.1) which is en- 
hanced in Na, N, and showing a correlation of Na with O. 
This group has the highest average barium of all the groups. 

High precision photometry was used to define 39 AGB 



Sbordone et al. (20111 calculated synthetic spectra 



based upon typical globular cluster abundances. These spec- 



stars in the |vL+07| sample. Their S(3839) and CH(4300) 
indices provided a good measure of the CN strength of 
the star. The majority of AGB stars were found to be 
CN-weak. The hypothesis that differing helium abundances 
could cause such an effect in mono-metallic clusters does 
not fit this result, as this AGB is associated with the most 
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metal-poor stars in the cluster and therefore will not have 
any helium enhancement. 
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